Fundamentals of Nano-Thermodynamics 



On the minimal lengtli scales, where temperature exists 

M. Hartmann, G. Mahler and 0. Hess 

Recent progress in the synthesis and processing of nano-structurcd materials and systems calls 
for an improved understanding of thermal properties on small length scales. In this context, the 
question whether thermodynamics and, in particular, the concept of temperature can apply on 
the nanoscale is of central interest. Here we consider a quantum system consisting of a regular 
chain of elementary subsystems with nearest neighbour interactions and assume that the total 
system is in a canonical state with temperature T. We analyse, under what condition the state 
factors into a product of canonical density matrices with respect to groups of n subsystems each, 
and when these groups have the same temperature T. In quantum systems the minimal group 
size rimin depends on the temperature T, contrary to the classical case. As examples, we apply 
our analysis to a harmonic chain and different types of Ising spin chains. For the harmonic chain, 
which successfully describes thermal properties of insulating solids, our approach gives a first 
quantitative estimate of the minimal length scale on which temperature can exist: This length 
scale is found to be constant for temperatures above the Debye temperature and proportional 
to below. We finally apply the harmonic chain model to various materials of relevance 

for technical applications and discuss the results. These show that, indeed, high temperatures 
can exist quite locally, while low temperatures exist on larger scales only. This has striking 
consequences: In quasi 1-dimensional systems, like Carbon-Nanotubes, room temperatures (300 
Kelvin) exist on length scales of 1 ^m, while very low temperatures (10 Kelvin) can only exist on 
scales larger than 1 mm. 
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I. INTRODUCTION 

Thermodynamics is among the most successfuUy and extensively appUed theoretical concepts in physics. Notwith- 
standing, the various limits of its applicability are not fully understood 13 ITlh . 

Of particular interest is its microscopic limit. Down to which length scales can its standard concepts meaningfully 
be defined and employed? 
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Besides its general importance, this question has become increasingly relevant recently since amazing progress in 
the synthesis and processing of materials with structures on nanometer length scales has created a demand for better 
understan ding of thermal properties of nanoscale devices, individual nanostructures and nanostructured materials 
(P; "s^; '3^ 1401) . Experimental techniques have improved to such an extent that the measurement of thermodynamic 
quantities like temperature with a spatial resolution on the nanometer scale seems within reach 10: 31.). In 
particular, the study of thermal transport in nano scale devices has made greate advances from an experimental 
point of view, while the definition of local temperatures and thus of temperature profiles remains unclear (Q). Carbon 
nanotubes play an important role in nano technology and nano scale thermal transport. According to a recent proposal 
ijlOl) a carbon nanotube filled with Gallium could be used as a thermometer, which has a spatial extension of only 
lOnm. 

Physical properties of smaller and smaller system sizes have been the subject of several experimental studies in 
recent years. Some examples are: Investigations of the size dependence of the melting point of atomic clusters l|34l) . 
the size dependence of the surface ferromagnetism of Pd fine particles l|3^ . existence of large critical currents of 
superconductors at nanosizes and an exceptional surface energy of Ag nanoparticles |2^. Also the system size needed 
to allow for thermalisation has been investigated. Electrons traveling through a quantum wire only become thermally 
distributed if the wire has a certain minimum length (,31.) . 

To provide a basis for the interpretation of present day and future experiments in nanoscale physics and technol- 
ogy and to obtain a better understanding of the limits of thermodynamics, it is thus indispensable to clarify the 
applicability of thermodynamical concepts on small length scales starting from the most fundamental theory at hand, 
i. e. quantum mechanics. In this context, one question appears to be particularly important and interesting: Can 
temperature be meaningfully defined on nanometer length scales? 

The existence of thermodynamical quantities, i. e. the existence of the thermodynamical limit strongly depends 
on the correlations between the considered parts of a system. For short range interactions, the effective interaction 
between a region in space and its environment becomes less relevant compared to its internal energy with increasing 
region size. The volume of the region grows proportional to its diameter cubed, while its surface grows proportional 
to its diameter squared. This scaling behavior is used to show that correlations between a region and its environment 
become neglegible in the limit of infinite region size and that therefore the thermodynamical limit exists (9: 26: 33). 

To explore the minimal region size needed for the application of thermodynamical concepts, situations far away 
from the thermodynamical limit should be analysed. On the other hand, effective correlations between the considered 
parts need to be small enough (12). The scaling of interactions between parts of a system compared to the energy 
contained in the parts themselves thus sets a minimal length scale on which correlations are still small enough to 
permit the definition of local temperatures. It is the aim of this chapter to study this connection quantitatively. 
Typical setups where the minimal length scale, on which temperature exists, becomes relevant are measurements of 
a temperature profile with very high resolution. 

In recent years, some attempts to generalise the theory of thermodynamics such that it applies to small systems, 
have been made ifl^ IT5 ITsI : Is^ . These attempts take into account that the thermodynamic limit does not apply 
to such small systems. On the other hand, they fail to quantitatively analyse the role of the interactions between 
adjacent systems. 

To analyse the question whether temperature exists locally, one first has to define when temperature is said to exist. 
We adopt here the convention, that a local temperature exists if the considered part of the system is in a canonical 
state. The motivation of this convention is twofold: 

Following the standard arguments of statistical mechanics, the canonical state is the maximum entropy state of 
an interacting system and hence its equilibrium state. The entropy of the system S is thus a function of its internal 
energy U and temperature can be defined via 1/T = dS/dU. 

A more practical motivation is the following: Temperature is measured via expectation values of observables, 
their second moments etc. The canonical distribution is an exponentially decaying function of energy characterised 
by one single parameter, temperature. This implies that there is a one to one mapping between temperature and 
the expectation values of observables, by which temperature is measured. Temperature measurements via different 
observables thus yield the same result, contrary to distributions with several parameters. In large systems with a 
modular structure, the density of states is a strongly growing function of energy l)37j) . If the distribution were not 
exponentially decaying, the product of the denstity of states times the distribution would not have a pronounced peak 
and thus physical quantities like energy would not have "sharp" values. 

To determine when temperature exists locally, we thus have to analyse when an equilibrium state exists locally, that 
is when parts of a large system are in a thermal equilibrium state. The entire system need not be in an equilibrium 
state itself. However, local equilibrium is not believed to exist if the total system is far from equilibrium (24). It is 
not very well understood, when local equilibrium exists without a global one and only a few exact results on this topic 
are known. In our approach, we thus focus on systems which are even globally in an equilibrium state. Nevertheless, 
our results are expected to be more general, provided the global state is close to an equilibrium one lf25l) . 



3 




FIG. 1 Groups of n adjoining subsystems are formed. 

Based on the above argument and noting that a quantum description becomes imperative at nanoscopic scales, 
the following approach appears to be reasonable: Consider a large quantum system, brought into a thermal state via 
interaction with its environment, divide this system into subgroups and analyse for what subgroup-size the concept 
of temperature is still applicable. 

This approach, when applied to concrete models, yields quantitative estimates of minimal length scales on which 
temperature can exist. Hence, we apply our concept to two classes of models: 

Spin chains have recently been subject of extensive studies in condensed matter physics and quantum information 
theory. Thus correlations and possible local temperatures in spin chains are of interest, both from a theoretical and 
experimental point of view ( 22; 39). We thus study spin chains with respect to our present purpose. 

Harmonic lattice models are a standard tool for the description of thermal properties of solids. We therefore apply 
our theory to a harmonic chain model to obtain estimates that are expected to be relevant for real materials and 
might be tested by experiments. 

This chapter is organized as follows: In section m we present the general theoretical approach which derives two 
conditions on the effective group interactions and the global temperature. We discuss the relation to correlations 
and several scenaria that might occur. In the following two sections we apply the general consideration to two 
concrete models and derive estimates for the minimal subgroup size. Section IIIII considers various types of Ising 
spin chains, which have a finite energy spectrum and show similar features as fermionic systems (e.g. correlated 
electrons). Section llVI then deals with a harmonic chain, a model with an infinite energy spectrum. Harmonic chains 
or lattices successfully describe thermal properties of insulating solids. In the following section[V]we therefore discuss 
the results for some real materiales modeled by harmonic chains. As a last point, we discuss the measurability of 
local tempeartures in section lvTl In the concluding scction lVlTl we finally compare the results for the different models 
considered and indicate further interesting topics. 

II. GENERAL THEORY 

We start by defining the Hamiltonian of our chain in the form (|27l) . 



where the index i labels the elementary subsystems. Hi is the Hamiltonian of subsystem i and the interaction 

between subsystem i and i + 1. We assume periodic boundary conditions. 

We now form Nq groups of n subsystems each as illustrated in figure^ (index i ^ {p, — l)n + j; /i = 1, . . . , Ng\ j — 
1, . . . , n) and split this Hamiltonian into two parts. 




(1) 



H = Ho + I, 



(2) 



where Hq is the sum of the Hamiltonians of the isolated groups. 




with 



n 




(3) 



i=i 



and / contains the interaction terms of each group with its neighbour group. 




(4) 
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We label the eigenstates of the total Hamiltonian H and their energies with the greek indices {ip, tp) and eigenstates 
and energies of the group Hamiltonian Hq with latin indices (a, b), 

H \ip) = \ip) and Ho \a) = Ea \a). (5) 

Here, the states \a) are products of group eigenstates 

Ng 

i«) = n i«M>, (6) 

where {H^ - /^„,^n+i) |a^) = E^\a^). E^j, is the energy of one subgroup only and Ea = J2^=i ^ti- 



A. Thermal State in the Product Basis 

We assume that the total system is in a thermal state with a density matrix which reads 

in the eigenbasis of H. Here, Z is the partition sum and (3 — (ksT)^^ the inverse temperature with Boltzmann's 
constant ks and temperature T. Transforming the density matrix (O into the eigenbasis of Hq we obtain 

{a\p\a) = / Wa{E)^dE (8) 

for the diagonal elements in the new basis. Here, the sum over all states \a) has been replaced by an integral over 
the energy. Eq is the energy of the ground state and Ei the upper limit of the spectrum. For systems with an 
energy spectrum that does not have an upper bound, the limit Ei ^ oo should be taken. The density of conditional 
probabilities Wa{E) is given by 

{\ip):E<E^<E+AE} 

where AE is small and the sum runs over all states \(p) with eigenvalues E^p in the interval [E, E + AE). 

To compute the integral of equation © we need to know the distribution of the conditional probabilities Wa{E). 
Following the proof of the central limit theorem for mixing sequences (20) as a guideline, one can show that the 
characteristic function of H does not change if a small subset of the is neglected. Truncating iJ in a suitable 
way, the characteristic function of the remainder of H factorises. If the remainder of H then fulfills the conditions 
for Lyapunov's version of the central limit theorem 1^ and Wa{E) converges to a Gaussian normal distribution in the 
limit of infinite number of groups Nq [T^ . 

lim Wa{E) = -=i- exp {- ' 7 V (10) 



^o^oo - ' ^Aa '\ 2A2 
where the quantities ai^d are defined by 

Ea = {a\H\a) - {a\Ho\a) (11) 
^ (a|i/» - {a\H\a)\ (12) 

Note that has a classical counterpart while A^ is purely quantum mechanical. It appears because the commutator 
[H,Hq\ is nonzero, and the distribution WaiE) therefore has nonzero width. 

The rigorous proof of equation (|10|l is given in |l3 ) and based on the following two assumptions: The energy of 
each group as defined in equation Q is bounded, i. e. 

(xI^mIx) < C (13) 

for all normalised states \x) and some constant C, and 

{a\H^\a) - {a\H\a)^ > No C (14) 
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for some constant C > 0. 

Equation H1U|I does not only hold for linear chains, but is equally valid for lattice models of two or three dimensions. 

In scenarios, where the energy spectrum of each elementary subsystem has an upper limit, such as spins, condition 
is met a priori. For subsystems with an infinite energy spectrum, such as harmonic oscillators, we restrict our 
analysis to states where the energy of every group, including the interactions with its neighbours, is bounded. Thus, 
our considerations do not apply to product states |a), for which all the energy was located in only one group or only 
a small number of groups. The number of such states is vanishingly small compared to the number of all product 
states. 

If conditions 1131) and (|14|l are met, equation ((SJ can be computed for Nq ^ 1: 



(a|/5|a> =^ exp 



-l3Ea 



erfc 



V2Aa 



erfc 



Ei^Ea+ PAI 

V2Aa 



(15) 



where Ea — Ea + Sa and erfc(x) is the conjugate Gaussian error function. The second term only appears if the energy 
is bounded and the integration extends from the energy of the ground state Eq to the upper limit of the spectrum Ei . 

Note that Ea is a sum of Na terms and that fulfills equation H14|) . The arguments of the conjugate error 
functions thus grow proportional to \/ Nq or stronger. If these arguments divided by Nq are finite (different from 
zero), the asymptotic expansion of the error function (Q) may thus be used for Nq ^ 1: 



exp I 



erfc (a;) 



— for X 

sjTTX 

exp i—x"^) 
2 + — !- for X 



(16) 



Inserting this approximation into equation (|15|l and using Eq < Ea < Ei shows that the second conjugate error 
function, which contains the upper limit of the energy spectrum, can always be neglected compared to the first, which 
contains the ground state energy. 

The same type of arguments show that the normalisation of the Gaussian in equation I|1U|) is correct although the 
energy range does not extend over the entire real line (—00,00). 

Applying the asymptotic expansion (|16() . equation (|15|l can be taken to read 



(a|p|a) = - exp 
for (£;o ~Ea-ea + (3Al) / {y/2NB A,) < and 



-f3 E, 



en - 



!3Al 



(17) 



exp 



{a\p\a) 



-PE„ 



{Ea 



E,f 



2A2 



^27rZ 



E„-Ea-ea + (3Ai 
A„, 



(18) 



for (£^0 -Ea-ea + 13 Al) / {V2N^ Aa) > 0. 

The off diagonal elements {a\p\b) vanish for \Ea — Eb\ > Aa + Ab because the overlap of the two distributions of 
conditional probabilities becomes negligible. For \Ea — Ei,\ < Aa + Af,, the transformation involves an integral over 
frequencies and thus these terms are significantly smaller than the entries on the diagonal. 



B. Correlations 

We now analyse when the thermal state is correlated in the product basis. To this end, one usually studies the 
off diagonal elements (a|p|6) of the density matrix in the product basis. However, since, in our case, the off diagonal 
elements can be estimated to be much smaller than the diagonal elements (a|/9|a), some information about correlations 
may be obtained from the latter ones. For an uncorrelated state, they factorise. 



Ng 

{a\p\a) = n^A* 



(19) 
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where = TTjr{a\p\a) and Trjj is the trace over aU groups except group number /i. Taking the logarithm of both 
sides of equation l(T?^ one gets 

Ng 

ln((aHa)) = ^ln(p^) (20) 

As follows from equations (|17|l and (|18(l this can only be the case for Ea + Ea — Eq > /3A^. Furthermore, the terms in 
the exponential on the rhs of equation (|17|l must be a sum of terms, each of which depends on the state of one group 
only. A term depending on the state ja^) of group fi must not depend on the state \a^+i) of the neighbouring group. 
In the models we consider here, these conditions coincide with the conditions for the existence of local temperatures 
because Ea and are sums of terms which depend on two neighbouring groups each. 



C. Conditions for Local Thermal States 

We now test under what conditions the diagonal elements of the (local) reduced density-matrices are also canonically 

distributed with some local inverse temperature IS^^^ for each subgroup ^ — 1, . . . , Nq- Since the trace of a matrix is 
invariant under basis transformations, it is sufficient to verify that they show the correct energy dependence. If we 

assume periodic boundary conditions, all reduced density-matrices are equal (pII^^ = Aoc for all /i) and the products 
of their diagonal elements are of the form (a|p|a) oc exp (—l3\ocEa)- We thus have to verify that the logarithm of the 
right hand side of equations H17|l and (|18|l is a linear function of the energy Ea, 

In ((a|/5|a)) « -Aoc Sa + c, (21) 

where Aoc and c are constants. Note that equation H21(l does not imply that the occupation probability of an eigenstate 
\ip) with energy E^ and a product state with the same energy Ea ~ E^ are equal. Even if Aoc and /3 are equal with 
very good accuracy, but not exactly the same, occupation probabilities may differ by several orders of magnitude, 
provided that the energy range is large enough. 

We exclude negative temperatures (/? > 0). Equation (|21|l can only be true for 

^^ + ''^-^' >f3-^, (22) 



as can be seen from equations (|17ll and (|18|) . Furthermore Ea and have to be of the following form: 

-£a + f ^2 « C^Ea + C2 (23) 

where ci and C2 are constants. Note that Ea and A^ need not be functions of Ea and therefore in general cannot be 
expanded in a Taylor series. 

To ensure that the density matrix of each subgroup /i is approximately canonical, one needs to satisfy (|23|l for each 
subgroup /i seperately; 



{Al_,+Al) + i^Al « c,E, + c, (24) 



where 



Ng 



{a\Itj.n,iin+i\a) with = > (25) 



Ng 



Al ^ {a\nl\a) - {a\n,\ay with A^ = ^ A^ (26) 



^ {a\n,-inu+n,n,-i\a) ~2{a\n,-i\a){a\n,\a). (27) 

ly—j.L — l 

Temperature becomes intensive, if the constant ci vanishes, 

|ci|«l Aoc = A (28) 



E) ■ (^\pM 




E 



FIG. 2 The product of the density of states rj{E) times the occupation probabilities {Lp\p\Lp) forms a strongly pronounced peak 
a.t E ^E. 

\ni{a\p\a)) 




E 



FIG. 3 ln((a|p|a)) for p as in equation l|15|l (solid line) and a canonical density matrix p (dashed line) for a harmonic chain. 

If this was not the case, temperature would not be intensive, although it might exist locally. 

It is sufficient to satisfy conditions 122|) and l|24|l for an adequate energy range i?min ^ < E^^x only. 

For large systems with a modular structure, the density of states is typically a rapidly growing function of energy 
Is^)- If the total system is in a thermal state, occupation probabilities decay exponentially with energy. The product of 
these two functions is thus sharply peaked at the expectation value of the energy E of the total system E+Eq —Tr{Hp), 
with Eq being the ground state energy. The energy range thus needs to be centered at this peak and large enough. 
On the other hand it must not be larger than the range of values E^ can take on. Therefore a pertinent and "save" 
choice for i?min and i?max is 



En 
Em 



\E 



M J mill ' 



1 E 
a Ng 
E 



+ 



a 



Ng 



En.) 
Ng J 

En.) 
Ng J 



(29) 



where a ^ 1 and E will in general depend on the global temperature. In equation 1)29(1 . [i?^] . and [E'^] denote 
the minimal and maximal values E^ can take on. 

Figure 121 shows the logarithm of equation (|15|l and the logarithm of a canonical distribution with the same /3 for a 
harmonic chain. The actual density matrix is more mixed than the canonical one. In the interval between the two 
vertical fines, both criteria (|^ and are satisfied. For E < Ei (|^ is violated and (|^ for E > E2- To allow 
for a description by means of canonical density matrices, the group size needs to be chosen such that Ei < -Bmin and 

E2 > £^max- 

For a model obeying equations H13() and (|14|l , the two conditions H22|) and 1(24(1 , which constitute the general result 
of this article, must both be satisfied. These fundamental criteria will now be applied to some concrete examples. 
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III. ISING SPIN CHAIN IN A TRANSVERSE FIELD 



In this section we consider an Ising spin chain in a transverse field. For this model the Hamiltonian reads 



H.. = -Bal 
2 



(30) 



where erf, erf and erf are the Pauli matrices. B is the magnetic field and and Jy are two coupling parameters. We 
will always assume B > Q. 

The entire chain with periodic boundary conditions may be diagonalized via succesive Jordan- Wigner, Fourier and 
Bogoliubov transformations (see appendix . The relevant energy scale is introduced via the thermal expectation 
value (without the ground state energy) 



E ■ 



uNg 

271 



Tfl ^ 



1 



where ujk is given in equation (jA9|) . The ground state energy Eq is given by 

uNg 



27r 



2 



(31) 



(32) 



Since Ng S> 1, the sums over all modes have been replaced by integrals. 

If one partitions the chain into Nq groups of n subsystems each, the groups may also be diagonalized via a Jordan- 
Wigner and a Fourier transformation (see appendix . Using the abbreviations 



K 



Jx H~ Jy 

2B 



and L 



Jx Jy 

2B 



the energy Ea reads 



Ng / 1\ 

= 25 ^ 5] [1 - if cos(fc)] [nl{^^)--] 



(33) 



(34) 



where k = irl / {n + 1) (I = 1, 2, . . . , n) and n'j:{fi) is the fermionic occupation number of mode k of group /i in the state 
I a). It can take on the values and 1. 

For the Ising model at hand one has Ea = for all states '"^ -^-^^i^ +1^^ a2 



while the squared variance reads 



A 



Ng 
11=1 



(35) 



with 



B' 



2 



L2 
2 



2 {K"^ - L^) 



-^sin2(fc) (nl{^,)-\ 



(36) 



-l_^3in^(,) (n> + l)4)]) 



where the n'j,(^) are the same fermionic occupation numbers as in equation H34|) . 

The conditions for the central limit theorem are met for the Ising chain apart from two exceptions. Condition l|13() 
is always fulfilled as the Hamiltonian of a single spin has finite dimension. As follows from equation H36(l , condition 
is satisfied except for one single state in the case where Jx — Jy [L = 0) and Jx = — Jy [K = 0) respectively. 
Therefore, for Ng ^ 1 the fraction of states where our theory does not apply is neglegible. 

We now turn to analyse conditions H22|l and H24|) . Since the spectrum of the Ising chain is limited, there is no 
approximation analog to the Debye approximation for the harmonic chain and A^^ cannot be expressed in terms of 



E, 



and E^. 



We therefore approximate H22|l and (|24|l with simpler expressions. 
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Let us first analyse condition 122(1 . Since it cannot be checked for every state \a) we use the stronger condition 

instead. This imhes that H22|l holds for all states \a). We require that H37|) is true for all states with energies in the 
range (P^ . It is hardest to satisfy for = E^in, we thus get the condition on n: 

^min ^0 

where e^in = E^^in/n and cq = Eo/inNc)- 

We now turn to analyse condition H24|) . Equation H3()|) shows that the do not contain terms which are proportional 

to Efj_. One thus has to determine, when the are approximately constant. This is the case if 

[A^l - [A^l 

g I MJmax L min l _ r^; l /ggN 

where [a;] max and [x]min denote the maximal and minimal value x takes on in all states \a). As a direct consequence, 
we get 

|ci| « 1 (40) 
which means that temperature is intensive. Defining the quantity = E^/n, we can rewrite (|39|) as a condition on 



n. 



2" 



max [^MJ min ( AA \ 

2 5 [eJ -[eJ . ^ ' 

L MJmax L MJmin 

where the accuracy parameter (5 ^ 1 is equal to the ratio of the Ihs and the rhs of H39|) . 

Since equation H39|l does not take into account the energy range (|29|) . its application needs some further discussion. 

If the occupation number of one mode of a group is changed, say from n'^ip) = to nf^{^j) = 1, the corresponding 
A2 differ at most by 45^ \k'^ _ 2,21 / _^ Qn the other hand, [A^l - [A^l . = 1^2 _ 2^21 r^-^^ ^^^^^ 

M *^ I I ' ^ ^ ' L /^J max L /^J min I I 

with the maximal A^ and the state with the minimal A^ thus differ in nearly all occupation numbers and therefore 
their difference in energy is close to \Eu\ — \En\ • . On the other hand, states with similar energies E,. also have 

'^■^ L MJmax L MJmm ' o 

a similar A^. Hence the A^ only change quasi continuously with energy and equation H39|l ensures that the A^ are 
approximately constant even on only a part of the possible energy range. 

To illustrate the scenaria that can occur, we are now going to discuss three special coupling models which represent 
extremal cases of the possible couplings in our model ((ISil . 



A. A coupling with constant Aai Jy ~ Q 

If one of the coupling parameters vanishes ( = Q or Jy = Q) , K — L and A^ = is constant. In this case 

only criterion H22|) has to be satisfied, which then coincides with (|38|l . 

Inserting expressions H31|) . (|32|) and ^A<6\ with = J and Jy = into condition H38|) . one can now calculate 
the minimal number of spins per group. Figure 0] shows rimin for weak coupling K = L = 0.1 and strong coupling 
K — L — 10 with a = 10 as a function of T/B. We choose units where Boltzmann's constant ks is 1. 

Note that, since A^ = const, condition (|38|1 coincides with criterion (|22l) . so that using (|38|) does not involve any 
further approximations. Figure [S] shows rinun according to eq. (|38l) as a function of the coupling parameter K and 
T/B for a — 10. The range with K < 1 represents the weak coupling regime, since the local level splitting is larger 
than the coupling strength. As one would expect, rimin grows with increasing coupling K. At low temperatures, rimin 
does not smoothly approach 1 for ^ 0. This unexpected feature originates from the Nq — > 00 limit, we consider 
here. 



B. Fully anisotropic coupling: = —J, 



If both coupling parameters are nonzero, the A^ are not constant. As an example, we consider here the case of 
fully anisotropic coupling, where Jx = — Jy, i. e. K = 0. Now criteria H38|) and (|41|l have to be met. 
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10"* 10"^ 10~* 10~^ 10^ 

FIG. 4 Log-log-plot of riniin according to eq. (I38I I for K = L — 0.1 (dashed line) and for K = L — 10 (solid line) as a function 
of T/B. a = 10 is defined in equation l)29|l . 




FIG. 5 nmin according to eq. 1381 1 as a function of the coupling parameter K and T/B. a = 10 is defined in equation 12911 . 
For K = 0, one has [A^l = L^, \Al] . = and [e^l = - \eJ . = B. 

' L MJ max L MJ min ^ M-lmax L MJinm 

We insert these results into (|41ll as well as (|31|l and l|32|l into (|38|l and calculate the minimal number of spins per 
group. Figure El shows n-aan according to criterion (|38l) and according to criterion (|41|) sepcratcly, for weak coupling 
L = 0.1 and strong coupling L = 10 with a = 10 and 5 — 0.01 as a function of T/B. 

In the present case, all occupation numbers n^.^^) are zero in the ground state. In this state, is maximal 
(A^ = B^ L^) as can be seen from (|36|) . Therefore criterion H38|) is equivalent to criterion (|22|l for K — 0. 

Figure [71 shows n^nin according to eq. (|38|l as a function of the coupling parameter K and T/B for a = 10. Again, 
L < 1 represents the weak coupling regime, where the local level splitting is larger than the coupling strength. As in 
the previous case, nmin is a monotonic increasing function of L but shows a non-smooth transition to the uncoupled 
(L = 0) case. The latter is again due to Nq oo limit. 

C. Isotropic coupling: Jx — Jy 

As a second example of models where both coupling parameters are nonzero, we consider the isotropic coupling 
case, where Jx = Jy, i. e. L = 0. Again, both criteria (|38|) and H41fl have to be met. 

The values of [A^]^^^^, [A^]^.^, [e^]^^^ and [e^],^;^ are given in equations dUIJ, dUSl) and jA13|). 

For the present model with L = and \K\ < 1 all occupation numbers n'^{^j) are zero in the ground state and 
thus A^ = 0. As a consequence, condition H38I) cannot be used instead of 1)22(1 . We therefore argue as follows: In 
the ground state — Eq/Nq = as well as A^ — and all occupation numbers n'j^ip) are zero. If one occupation 
number is then changed from to 1, A^ changes at most by 4 B^ K"^ /{n + 1) and E^ changes at least by 2 _B (1 — jiirl). 
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FIG. 6 Log-log-plot of riniin according to eq. 13811 for L = 0.1 (dashed line) and for L = 10 (solid line) and rimin according to 
eq. H41^ for L = 0.1 (dash - dot line) and for i = 10 (gray line) as a function of T/B. A" = 0, a = 10 and 5 — 0.01. a and 5 
are defined in equations 1)29^ and H41|l respectively. 




FIG. 7 riniin according to eq. 13811 and eq. 14111 as a function of the coupling parameter L and T/B. q = 10 is defined in 
equation 129|l . 



Therefore (|22(l will hold for all states except the ground state if 

">2i3/3Y^ (42) 

If \K\ > 1, occupation numbers of modes with cos(fc) < l/|ii^| are zero in the ground state and occupation numbers 
of modes with cos(fc) > are one. for the ground state then is [A^]^^^ « [A^]^^^/2 and is a good 

approximation of condition H22I) . 

Inserting these results into iglj as well as (j2U and (j221) into (jSHl for \K\ > 1 and using |g2l for \K\ < 1, the 
minimal number of spins per group can be calculated. Figure |S1 shows Uynin according to criteria (|42l) and H41|) for 
weak coupling if = 0.1 and according to criteria H38() and (|41|l for strong coupling K — 10 with a = 10 and S = 0.01 
as a function of T/B. 



IV. HARMONIC CHAIN 

As a representative for the class of systems with an infinite energy spectrum, we consider a harmonic chain of Nq ■ n 
particles of mass m and spring constant ^JmLOQ. In this case, the Hamiltonian reads 

= ^pI + y^oQ' (43) 
= -mujlq.iqi+i, (44) 
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FIG. 8 Log-log-plot of nmin according to eq. (I42I I for if = 0.1 (dashed line), accor ding to eq. igHJ for if = 10 (solid line) 
and according to eq. H41|l for K — 0.1 (dash - dot line) and for if = 10 (gray line) as a function of T/B. L = 0, a = 10 and 
6 = 0.01. a and 5 are defined in equations 1)29^ and l)41|l respectively. 



where pi is the momentum of particle i and qi the displacement from its equilibrium position i ■ oq with oq being the 
distance between neighbouring particles at equilibrium. We devide the chain into Nq groups of n particles each and 
thus get a partition of the type considered in section ITll 

The Hamiltonian of one group is diagonalised by a Fourier transform and the definition of creation and anihilation 
operators and Ofe for the Fourier modes (see appendix IB|| . 

Ng / -, X 



ti=l k 



where k = nl/{ao {n + 1)) (/ = 1,2, ... ,n) and the frequencies tJk are given by = Aojq sin^(fca/2). n'^{pi) is the 
occupation number of mode k of group in the state \a). We chose units, where h= \. 

To see that this model satisfies the condition (|14|) one needs to express the g roup interdiction \^{(l^n-} Qfin-\-i 

) in terms 

of and a^, which yields — for all /j, and therefore 

Ng 

^l = T.^l (46) 



where A^, the width of one group interaction, reads 




^cos2(^) {^p + ^ 



(47) 



A^ has a minimum value since all Uk > and all nip > 0. In equation (|47|l . k labels the modes of group /i with 
occupation numbers rik and p the modes of group /i + 1 with occupation numbers irip. The width A^ thus fulfills 
condition ((Til) . 

Since the spectrum of every single oscillator is infinite, condition H13|l can only be satisfied for states, for which the 
energy of the system is distributed among a relevant fraction of the groups, as discussed in section HTl 

The expectation values of the group interactions vanish, = 0, while the width s A^ depend on the occupation 
numbers Uk and therefore on the energies i?^. We thus apply the conditions (|22() and (|24(l . To analyse them, we make 
use of the continuum or Debye approximation l)23j) . requiring n 3> 1, ao I, where I — nao, and the length of the 
chain to be finite. In this case we have LUk = v k with the constant velocity of sound v ~ LUoaQ and cos(fc ao/2) « 1. 
The width of the group interaction H47f) thus translates into 

Al = ^E,,E,,+, (48) 



13 



where n + 1 « n has been used. The relevant energy scale is introduced by the thermal expectation value of the entire 
chain 

E = NGnkBei-^j (49) 
and the ground state energy is given by 

E = Nankse (|) ' | d. ^ ^^^^ (50) 

We first consider criterion H22I) . For a given Ea — J^fi^t^' the squared width is largest if all E^^ are equal, 
X = E^ Vfi. Thus H22(l is hardest to satisfy for that case, where it reduces to 

.-^-|.'>0^ ,51, 

Equation (|51|l sets a lower bound on n. For E > Eq, the bound is strongest for high energies, x = a{E / Nq) + {Eq / Nq) , 
while for £' < i?o, it is strongest for low energies, x = [E/aNc) + (Eq/Ng). Since condition H24fl is a stronger criterion 
than condition (|22|l at i? > i?o, we only consider 1)51(1 for low energies, where it reads 

Q a / e \^ , , 

n>4-- - + eo (52) 
1 e \a J 

with e = E/inNckBO) and eo = Eo/inNokBO) = 1/4. 

To test condition (I24|l we take the derivative with respect to E"^ on both sides, 

- ^E,., +E,^,-2—^+-—^c, (53) 

where we have seperated the energy dependent and the constant part on the Ihs. (|53|l is satisfied if the energy 
dependent part is much smaller than one, i.e. 

^(£;.-i+i?,+i-2^)<<5«l, (54) 

Taking E^^i and equal to the upper bound in equation 129|l . this yields 

" > — y e, (55) 

where the "accuracy" parameter (5 <C 1 quantifies the value of the energy dependent part in ^ 
Since the constant part of the Ihs of H53|l satisfies 

< — I — -— I (56) 



n 



temperature is intensive. 

Inserting equation 1(49(1 into equation 1(52(1 and ((55(1 one can now calculate the minimal n for given (5, a, and T. 
Figure El shows rimin for a = 10 and 5 — 0.01 given by criterion 152|) and 1(55(1 as a function of T/O. 
For high (low) temperatures n^m can thus be estimated by 

2 a/ 5 for T > 6 

(3a/27r2) {Q/Tf for T < 6 ^ ' 

Figure [Till shows rimin for a = 10 and i5 = 0.01 given by criterion ((52(1 and 1(55(1 as a function of the temperature T 
and the Debye temperature 0. As to be expected, nmin is a monotonic increasing function of 0. 
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FIG. 9 Log-log-plot of rimin according to eq. (I52II (dashed line) and rimin according to eq. I55I I (solid line) for a = 10 and 
5 — 0.01 as a function of T/O for a harmonic chain. S and a are defined in equations I55II and 1291 . respectively. 
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FIG. 10 Log-log-plot of Wniin according to eq. I|52|l and eq. I|55|l for a = 10 and 5 = 0.01 as a function of T and the Debye 
temperature for a harmonic chain. S and a are defined in equations 15511 and 1291 . respectively. 

V. ESTIMATES FOR REAL MATERIALS 

Thermal properties of insulating solids can successfully be described by harmonic lattice models. Probably the 
most well known example of such a successfull modeling is the correct prediction of the temperature dependence of 
the specific heat based on the Debye theory |2^. We therefor expect our approach for harmonic lattice models to 
give feasible estimates for real existing materials. 

The minimal length scale on which intensive temperatures exist in insulating solids, is thus given by 



(58) 



where ao is the lattice constant, the distance between neighbouring atoms. 
Let us consider some materials as examples: 



A. Silicon 

Silicon is used in many branches of technology. It has a Debye temperature of ~ 645 K and its lattice constant 
is ao ~ 2.4 A. Using these parameters, figure ITTl shows the minimal length-scale on which temperature can exist in a 
one-dimensional silicon wire as a function of global temperature. 
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B. Carbon 

Recently, carbon became increasingly important in the fabrication of nano-structured devices (0;0)- It has a Debye 
temperatm'e of « 2230 K and a lattice constant of oq ~ 1.5 A. 

Figure [T21 shows the minimal length-scale on which temperature can exist in a one-dimensional carbon device as a 
function of global temperature. 

Carbon Nanotubes have diameters of only a few nanometers. Figure 1121 thus provides a good estimate of the 
maximal accuracy, with which temperature profiles in such tubes can be meanigfuUy discussed (p). 



C. Gallium, Arsen and Indium 

As further representatives of materials which are of central relevance in todays optoelectronics and photonics, 
Figure El shows nmin for a = 10 and 6 — 0.01 given by criterion (|52|l and H55I) as a function of T/Q. we consider the 
semiconductors Gallium (9 w 320 K, oq ~ 2.4 A), Arsen (9 w 282 K, oq « 3.2 A) and Indium (9 w 108 K, ao w 3.3 A). 
Figure [T^ shows ^min for those materials. 



VI. LOCAL TEMPERATURE MEASUREMENT 

Having determined the limits on the existence of local temperatures, one is of course interested in whether and 
under what conditions temperature can be measured locally. We thus discuss this question in the present section. 

The standard setup for a temperature measurement consists of the system or a part of it, whose temperature is to be 
measured and the thermometer. The latter must be significantly smaller than the system and the interaction between 
the system and the thermometer must be weak, so that the exchange of energy between both cannot significantly 
alter the proper energy of the system. The thermometer furthermore should have the property, that there exists a 
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FIG. 13 /min as a function of temperature T for Gallium (dashed line), Arsen (solid line) and Indium (dash-dotted line) 



FIG. 14 Measurement setup 

one to one mapping between its temperature and another quantity, which is easily accessible to the observer, like the 
height of a meniscus or the pressure of a gas. 

If temperature exists for the system or its part of interest, the respective part is in a thermal equilibrium state. 
The thermometer then relaxes into an equilibrium state of the same temperature and thus measures the latter. From 
this point of view it is unclear whether the thermometer would measure the local or global equilibrium state. 

However, if one views the thermometer as a part of the system, the situation becomes rather obvious. If there 
exist an equilibrium state in the part of the system which contains the thermometer, this part can be described by a 
canonical density matrix. Since the interaction between the system (without the thermometer) and the thermometer 
is very weak this density matrix factorises into two canonical ones of the same temperature, one for the system and one 
for the thermometer. Therefore, the thermometer measures the temperature on a length scale where the "reduced" 
density matrix is closest to a canonical one. 

For a system in global equilibrium, the thermometer will thus always measure the global temperature. On the other 
hand, if equilibrium exists only locally, the thermometer measures the local temperature. 

The length scales derived in this chapter are thus only measurable, by standard methods, in scenarios of global 
non-equilibrium. It is therefore important to note that our results, although they are derived for global equilibrium 
scenarios, should also apply in more general settings. 

Whether local temperature can be measured in global equilibrium cases, e.g. via a measurement of occupation 
probabilities of the phonon spectrum, is an interesting subject of future research. 

VII. DISCUSSION AND CONCLUSIONS 

In summary, we have analysed the limits of the existence of temperature on the nanoscalc. To this end, we have 
considered a linear chain of particles interacting with their nearest neighbours. We have partitioned the chain into 
identical groups of n adjoining particles each. Taking the number of such groups to be very large and assuming the 
total system to be in a thermal state with temperature T we have found conditions (equations (|22|l and (|24|l 'l. which 
ensure that each group is approximately in a thermal state. Furthermore, we have determined when the isolated 
groups have the same temperature T, that is when temperature is intensive. 

The result shows that, in the quantum regime, these conditions depend on the temperature T, contrary to the 
classical case. The characteristics of the temperature dependence are determined by the energy dependence of Aq. 
The low temperature behavior, in particular, is related to the fact that Aq has a nonzero minimal value. This fact 
does not only appear in the harmonic chain or spin chains but is a general quantum feature of most systems of that 



part of the system to be measured 
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topology. The commutator [H, Hq] is nonzero and the ground state of the total system is energetically lower than the 
lowest product state, therefore aa is nonzero, even at zero temperature l|a l2lt ISOt l39|) . 

We have then applied the general method to a harmonic chain and several types of Ising spin chains. For concrete 
models, the conditions H22I) and 124|) determine a minimal group size and thus a minimal length scale on which 
temperature may be defined according to the temperature concept we adopt. Grains of size below this length scale 
are no more in a thermal state. Thus temperature measurements with a higher resolution should no longer be 
interpreted in a standard way. 

The most striking difference between the spin chains and the harmonic chain is that the energy spectrum of the 
spin chains is limited, while it is infinite for the harmonic chain. For spins at very high global temperatures, the total 
density matrix is then almost completely mixed, i. e. proportional to the identity matrix, and thus does not change 
under basis transformations. There are thus global temperatures which are high enough, so that local temperatures 
exist even for single spins. 

For the harmonic chain, this feature does not appear, since the size of the relevant energy range increases indefinitely 
with growing global temperature, leading to the constant minimal length scale in the high energy range. 

Our approach should be seen as a first estimate of these minimal length scales. It is thoroughly derived, nonetheless, 
the treatement may still be made more precise in some points. On the one hand, an exact treatement of the off- 
diagonal elements of the density matrix would give more justification to the results. On the other hand, the choice of 
the relevant energy range H29I) could be made more precise, taking into account the shape of the peak (see figure |2J|. 

The results we have obtained for the harmonic chain model allow to make predictions for real existing materials. 
We have thus considered some materials, which play an important role in fields of current technological interest. 

We have then discussed the measurability of local temperatures in a standard scenario and found that a thermometer 
always measures temperature on a length scale, where the "reduced" density matrix is closest to a canonical one, i.e. 
closest to an equilibrium state. Hence a thermometer measures the global temperature as long as the total system is 
in equilibrium, while it measures the local one if equilibrium exist only locally. 

The length scales, calculated in this paper, should also constrain the way one can meaningfully define temperature 
profiles in non-equilibrium scenarios |28l). Although, our approach only considers cases where even the total system 
is in an equilibrium state, the results are expected to be more general. 

A typical non-equilibrium scenario would be a heat conduction experiment, where the two ends of a wire are 
conected to two heat baths at different temperatures. In this setup, a temperature gradient would build up in the 
wire. If the gradient is small, the density matrix of the wire would be somewhat similar to a canonical one. Applying 
our approach to this quasi-equilibrium case, would then yield results of the same order of magnitude as the present 
treatment of the full equilibrium case. 



APPENDIX A: Diagonalization of the Ising Chain 

The Hamiltonian of the Ising chain is diagonalised via Jordan- Wigner transformation which maps it to a fermionic 
system. 



c, = (7 



.Z \ I ' z 



\j<i 



The operators Ci and c\ fulfill fermionic anti-commutation relations 

{c„Cj} = 



and the Hamiltonian reads 



H = B 



c, -1] -K 



' E + h.c.) - L ^ (44^, + h.c. 



(Al) 



(A2) 



(A3) 



with K = {J^ + Jy)/{2B) and L — {J^ — Jy)/(2B). In the case of periodic boundary conditions a boundary term 
is neglected in equation ljA3|l . For long chains {nNc oo) this term is supressed by a factor {nNc)~^- The 
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Hamiltonian now describes Fermions which interact with their nearest neighbours. As for the bosonic system, a 
Fourier transformations maps the system to noninteracting fermions. For the whole chain with periodic boundary 
conditions 



1 



Ci ) \/nNG 



with k = {2t:1) I [uNg) where ^ = 0, ±1, . . . , ±{nNG - 2)/2, {nNG)/2 for uNg even, and 



'3 '> = 



(A4) 



(A5) 



with k = (ttI) /{n + 1) and (Z = 1, 2, . . . , n) for one single group. 

In the case of periodic boundary conditions, fermion interactions of the form djjd^j, and dfed-fe remain. Therefore, 
one still has to apply a Bogoliubov transformation to diagonalize the system, i.e. 

dk = Ukbk + ivkh^_i, (A6) 
— 1. With the definitions Uk = cos(8a;/2) and Vk = sin(8A;/2) the interaction 



where Uk = u-k, Vk : 
terms disappear for 



-v-k and 



cos(efc) 



I — K cos k 



^[1 - A"cosA:]2 + [LsinA:]2 



(A7) 



In the case of the finite chain of one group, the Bogoliubov transformation is not needed, since the corresponding 



terms are of the form dj^dj. and dkdk and vanish by virtue of equation (|A2I 
The Hamiltonians in the diagonal form read 



H = J2^k[blbk 



where the frequencies are 



LUk ^ 2Bv/[l -Kcosk]^ + [Lsin/c]2 
with k — {2ttI) / (uNg) for the periodic chain and 

LUk = 2B(1 - Kcosk) 



(AS) 



(A9) 



(AlO) 



with k ~ {Til)/{n + 1) for the finite chain. 

For the finite chain the occupation number operators may also be chosen such that LOk is always positive. Here, the 
convention at hand is more convenient, since the same occupation numbers also appear in the group interaction and 
thus in A,,. 



1. Maxima and minima of iSf, and 

The maximal and minimal values of are given by 



MJmax L — 



for \K\ < 1 and by 



max 

' 1 

/^Jmin 



nB, 



\/ — 1 + arcsin 



\K\ 



(All) 



(A12) 



for I if I > 1, where the sum over all modes k has been approximated with an integral. 



The maximal and minimal values of A are given by 



L A^J mm 



max (if 2 L2) 
min {K^,L^) 



(A13) 
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APPENDIX B: Diagonalization of the Harmonic Chain 

The Hamiltonian of a harmonic chain is diagonaUsed by a Fourier transformation and the definition of creation and 
anihilation operators. 

For the entire chain with periodic boundary conditions, the Fourier transformation reads 

1 ^ 1 f Mfe exp(mofcj) 1 

with k = 2nl/{aonNG) and {I = 0,±1,..., 

±{nNG — 2)/2, {nNG)/2, where hNq has been assumed to be even. 

For the diagonahsation of one single group, the Fourier transformation is 

{2}-y^i:{:}^i..(«o«) (B2) 

with k = nl/{aQ (n + 1)) and (/ = 1,2, . . . ,n). 

The definition of the creation and anihilation operators is in both cases 



41 1 



muJkUk { , \ivk] (B3) 



I I ""^K'^K \ I 

flfc J ^/2rwjJk \ t + 

where the corresponding Uk and Vk have to be inserted. The frequencies ujk are given by = Awq s\v?{kao/2) in both 
cases.. 

The operators a\ and satisfy bosonic commutation relations 

[a/c,ap] = 

[ afc,aj, ] = 6kp (B4) 



and the diagonahsed Hamiltonian reads 



H = Y,^k{alak + ^ (B5) 



References 

[1] M. Abramowitz and I. Stegun. Handbook of Mathematical Functions. Dover PubL, New York, 9th edition, 1970. 

[2] A. E Allahverdyan and Th. M Nieuwenhuizen. Extraction of Work from a Single Thermal Bath in the Quantum Regime. 

Phys. Rev. Lett, 85:1799-1802, 2000. 
[3] A. E. Allahverdyan and Th. M. Nieuwenhuizen. Testing the Violation of the Clausius Inequality in Nanoscale Electric 

Circuits. Phys. Rev. B, 66:115309, 2002. 
[4] J. Aumentado, J. Eom, V. Chandrasekhar, P.M. Baldo, and L.E. Rehn. Proximity effect thermometer for local electron 

temperature measurements on mesoscopic samples. Appl. Phys. Lett., 75:3554, 1999. 
[5] P. Billingsley. Probability and Measure. John Wiley & Sons, New York, 3rd edition, 1995. 

[6] D. Cahill, W. Ford, K. Goodson, G. Mahan, A. Majumdar, H. Maris, R. Merlin, and S. Phillpot. Nanoscale thermal 

transport. J. Appl. Phys., 93:793, 2003. 
[7] M.S. Drcssclhaus, editor. Carbon Nanotubes, volume 80 of Topics in Applied Physics. Springer, Berlin, 2001. 
[8] M.S. Drcsselhaus, G. Dresselhaus, and P.C. Eklund. Science of Fullerenes and Carbon Nanotubes. Academic Press, San 

Diego, 1996. 

[9] M.E. Fisher. The Free Energy of a Macroscopic System. Arch. Rati. Mech. Anal., 17:377, 1964. 

[10] Y. Gao and Y. Bando. Carbon nanothcrmonieter containing gallium. Nature, 415:599, 2002. 

[11] J. Gcmmcr, A. Ottc, and G. Mahler. Quantum approach to a derivation of the second law of thermodynamics. Phys. Rev. 
Lett, 86:1927-1930, 2001. 

[12] M. Hartmann, J. Gemmer, G. Mahler, and O. Hess. Scaling behavior of interactions in a modular quantum system and 

the existence of local temperature. Euro. Phys. Lett, 65:613-619, 2004. 
[13] M. Hartmann, G. Mahler, and O. Hess. Existence of Temperature on the Nanoscale. Phys. Rev. Lett, (accepted), 2004. 

quant-ph/0312214. 

[14] M. Hartmann, G. Mahler, and O. Hess. Gaussian Quantum Fluctuations in Interacting Many Particle Systems, math- 
ph/0312045, 2004. acc. for publ. in Lett. Math. Phys. 



20 



[15] M. Hartmann, G. Mahler, and O. Hoss. Local Versus Global Thermal States: Correlations and the Existence of Local 

Temperatures, quant-ph/0404 1 64 , 2004. 
[16] M. Hartrnarui, G. Mahler, and O. Hess. Spectral densities and partition functions of modular quantum systems as derived 

from a central limit theorem, cond-mat/0406100, 2004. 
[17] T.L. Hill. A Different Approach to Nanothermodynamics. Nano Lett, 1(5):273, 2001. 

[18] T.L. Hill. Extension of Nanothermodynamics to Include a One-Dimensional Surface Excess. Nano Lett., 1(3):159, 2001. 
[19] T.L. Hill. Perspective: Nanothermodynamics. Nano Lett, 1(3):111, 2001. 

[20] LA. Ibargimov and Y.V. Linnik. Independent and Stationary Sequences of Random Variables. Wolters-Noordhoff, Gronin- 

gen/Nethcrlands, 1971. 

[21] A. N. Jordan and M. Biittiker. Entanglement Energetics at Zero Temperature. Phys. Rev. Lett, 92:247901, 2004. cond- 
mat/0311647. 

[22] M. Kenzelmann, R. Coldea, D.A. Tennant, D. Visser, M. Hofmann, R. Smeibidl, and Z. Tylczynski. Order-to-disorder 
transition in the XY-like quantum magnet Cs2CoC14 induced by noncommuting applied fields. Phys. Rev. B, 65:144432, 
2002. 

[23] Ch. Kittel. Quantum Theory of Solids. Wiley, New York, 1963. 

[24] H. J. Kreuzer. Nonequilibriurn Thermodynamics and its statistical foundation. Clarandon Press, Oxford, 1981. 

[25] R. Kubo, M. Toda, and N. Hashitsume. Statistical Physics II: Nonequilibriurn Statistical Mechanics. Number 31 in 

Solid-State Sciences. Springer, Berlin, Heidelberg, New- York, 2. edition, 1991. 
[26] J.L. Lebowitz and E.H. Lieb. Existence of Thermodynamics for Real Matter with Coulomb Forces. Phys. Rev. Lett., 

22:631, 1969. 

[27] G. Mahler and V.A. Weberrufi. Quantum Networks. Springer, Berlin, Heidelberg, 2. edition, 1998. 

[28] M. Michel, M. Hartmann, J. Gernrner, and G. Mahler. Fourier's Law confirmed for a class of small quantum systems. 

Euro. Phys. J. B, 34:325-330, 2003. 
[29] K.K. Nanda, A. Maisels, F.E. Kruis, H. Fissan, and S. Stappert. Higher Surface Energy of Free Nanoparticles. Phys. Rev. 

Lett, 91:106102, 2003. 

[30] Th. M. Nieuwenhuizen and A. E. AUahverdyan. Statistical Thermodynamics of Quantum Brownian Motion: Construction 

of Perpetuurn Mobile of the Second Kind. Phys. Rev. E, page 036102, 2002. 
[31] H. Pothier, S. Gueron, N.O. Brige, D. Estevc, and M.H. Devoret. Energy Distribution Function of Quasiparticles in 

Mesoscopic Wires. Phys. Rev. Lett, 79:3490, 1997. 
[32] A.K. Rajagopal, C.S. Pande, and Sumiyoshi Abe. Nanothermodynamics - A generic approach to material properties at 

the nanoscale. cond-mat/0403738, 2004. 
[33] D. Ruelle. Statistical Mechanics. W.A. Benjamin Inc., New York, 1969. 

[34] M. Schmidt, R. Kusclie, B. von Issendorf, and H. Haberland. Irregular variations in the melting point of size-selected 

atomic clusters. Nature, 393:238, 1998. 
[35] K. Schwab, E.A. Henriksen, J.M. Worlock, and M.L. Roukes. MeEisurement of the quantum of thermal conductance. 

Nature, 404:974-976, 2000. 

[36] T. Shinohara, T. Sato, and T. Taniyama. Surface Ferromagnetism of Pd Fine Particles. Phys. Rev. Lett, 91:197201, 2003. 
[37] R.C. Tolman. The Principles of Statistical Mechanics. Oxford Univ. Press, London, 1967. 

[38] J. Varesi and A. Majumdar. Scanning Joule expansion microscopy at nanometer scales. Appl. Phys. Lett., 72:37, 1998. 
[39] X. Wang. Threshold temperature for pairwise and many-particle thermal entanglement in the isotropic heisenberg model. 

Phys. Rev. A, 66:044305, 2002. 
[40] C.C. Williams and H.K. Wickramasinghe. Scanning thermal profiler. Appl. Phys. Lett., 49:1587, 1986. 



